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We present an approximate global solution of Einstein's equations which des- 
cribes a simple star model: a self-gravitating perfect fluid with constant density 
and rigid rotation (Born). In order to do this we carry out a second order post- 
minkowskian approximation of the problem, and then we use a slow rotation 
subapproximation. The resulting metric depends on three arbitrary constants: 
density, rotational velocity and stellar radius in the non-rotating limit. Mass, 
angular momentum, quadrupole and other moments are expressed in terms of 
these three constants. As a side result we show that this type of fluid cannot 
be a source of the Kerr metric. 



1 Introduction 

In this work we present an approximate solution of the Einstein equations 
describing a global model for the gravitational field generated by a bounded, self- 
gravitating stationary and axisymmetric body rotating rigidly with constant angular 
velocity. The exterior metric corresponds to an asymptotically flat stationary, ax- 
isymmetric vacuum solution, while the interior is a perfect-fluid solution with the 
above properties. Both solutions (exterior and interior) are matched at the surface of 
zero pressure using the Lichnerowicz matching conditions, i. e., the global solution is 
C 1 at the boundary. To carry out the matching we use global harmonic coordinates, 
since these coordinates provide the necessary arbitrary constants. 

1.1 Generic expresions of the metric 

Let us consider a stationary metric with axial symmetry. As usual, the metric 
tensor will be referred to a system of spherical coordinates {t, r, 9, <p} adapted to the 
symmetries and to the Papapetrou structure, so that the radial and polar coordi- 
nates are defined on the 2-D surfaces orthogonal to the surfaces of transitivity. For 
convenience, we shall use an orthonormal cobasis, with respect to the Minkowskian 
metric, associated to these coordinates: 
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g = g 0Q LJ°(g)LJ + g^^fgiiv 3 +lj 3 <E)iv°) + g 33 uj 3 (g)LJ 3 

+ g n u) 1 ®^ 1 +3i2 (uJ 1 ®^ 2 +uj 2 ®u> 1 ) + g 22 u 2 ®ijj 2 (1) 

uj°=dt , u) 1 = dr , uj 2 =rdd , u> 3 = rsin8d(p (2) 

However, for computational reasons, we shall write the metric tensor g with respect 
to another tensor basis 

g = 5oot 00 + 5o 3 t 03 + f d t d + jot + fit 1 + / 2 t 2 (3) 

t d = Sn dx l ®dx j 



t 00 = u/W° I t° = (Sij - 3e i e j )dx i <2>dx j 

t 03 = u°<E) u) 3 + u 3 <E) lj° | t 1 = (he-j + kjdjdx 1 ® dx° 

? = (kikj — mjmijdx' 1 ® dx^ 



(4) 



where {x 1 } are the "cartesian" coordinates associated to the spherical ones {r, 6, ip} 
and the covectors {fcj,ej,mj} arc defined as 

d(r sin 8) — h dx % , dx 3 = e% dx % , lu 3 = r sin 8 dip = rm dx l (5) 

that is to say, they represent respectively the unitary radial, axial and azimuthal 
cylindrical covectors. The tensor basis {t 00 , t 03 , t d , t°, t 1 , t 2 } appears quite naturally 
in the linear and higher order approximations of Einstein's equations, which require 
a Poisson-likc equation to be solved, which in turn requires the introduction of the 
tensor spherical harmonics (we shall see that they are directly related to the tensor 
basis we are considering). Notice that the tensor t d is pure trace, while {t ^ 1 ^ 2 } 
are traceless (with respect to the flat 3-dimcnsional metric). 

The relationship between the metric components on the spherical co-basis (2) 
and the coefficients of the metric with respect to the "pure trace-traceless" tensor 
basis (4) is: 

fd= ^(.9ll + .922 +333) 

fo = I [911 (1 - 3 cos 2 8) + 3.912 sin 28 + .g 22 (l - 3 sin 2 8) + g 33 ] 

? (6) 

h = 2^ gn ~ - 922 ) sin 20 + 912 cos 20 

h = 7; (.911 sm2 + .912 sin 28 + g 2 2 cos 2 8 - g 33 ) 
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1.2 Post— minkowskian approximation 

In this section we shall briefly describe the post-minkowskian approximation 
using the harmonic condition (we could also use the so called q-armonic condition 
[1]). We begin by defining the deviation of the metric in "cartesian" coordinates 

h a /3 ■= g a p - Va0 (7) 

Let us consider the Einstein equations and harmonicity conditions 

R a p = 8ir (T a p - \g a pT) = 8irT a/3 

(8) 

(we assume G = c = 1). As it is known the Ricci tensor and the l.h.s. of the harmonic 
condition can be written in terms of the deviation h a p as follows: 



R 



ap 



d p h 



11 pa 



+ N a p (9) 



d p h 



-jjlpa h 



Hn 



(10) 



where N a p and H a are at least quadratic in h a p. We also use the notation h = 
"f] a ^h a p , d p = n px d\ and □ = n Xp d\^ (usual d'Alembertien operator). Consequently 
the Einstein equations using harmonic coordinates can be written as: 



□ Kp = -l6vT a p + 2 [N a p - d {a Hp } ] 

2 ^P a ^) " 



(ii) 



d p {h pa 

and then the equations for the linear stationary approximation can be written as 
follows (where A is the flat laplacian operator): 



Ah 



Qk 



(1) 

ap 



-16wT 



(i) 



a/3 



l ka 



Vka 







(12) 



where, as shown below, the tensor associated to the energy-momentum tensor 
requires only the flat Minkowski metric. Once a solution of these equations is fixed 
we can solve the postlinear approximation from the equations 



LAn ap — lV>1T1 ap 



2N 



(2) 
ap 



d h (2) 



(2) 



,(2) 
L ka 



Vka 



(13) 
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(2) (2) (2) (1) 

where , and H a ' are constructed from the linear deviation h a p. Note 

that using the expressions (9,10) it obviously follows that and H a 2 ^ may be 

obtained in an algorithmic way easily implemented on a computer. One only needs 
evaluate this expressions for the metric gj^ = r/ a fj + h^p 



2 Exterior solution 

2.1 Linear approximation (exterior solution) 

The general asymptotically flat solution of equations (12) in vacuum (exterior) 
with the assumed symmetries is 



n—0 n—1 n=0 

+ E ^r E « + E (1 4 ) 



n=l 

where the notation: 



En = -|nD„ + ^nH° + H* 

F„ = l -n{n - 1)H° + (n - l)Hj, - ^H* 



has been used, and 



(15) 



f T„ = P„(cos0)t 00 f D„ = P„(cos0)t d , H° = P n (cos0)t° 
1 Z„ ^P r i(cos0)t 03 \ Hi sP^cos^t 1 , = p2(cos0)t a 

arc the axially symmetric tensor spherical harmonics. The constants Mn^ and Jn^ 
represent the linear intrinsic multipole moments of Geroch and Hansen [4], while 

the constants ^ and are related to the stress moments of the source in this 
linear approximation [6]. The constants A and B could be eliminated by choosing 
an appropriate gauge, but they are necessary for the matching problem with the 
Lichnerowicz conditions [8]. 

In the following wc ncg lcct the multipolcs M^ ] , A {1) , B^ ] and following. 
We also assume equatorial symmetry which eliminates the odd (M, A, B) and even 
J multipole moments. This is coherent with the approximations done in the last 
Section. 



An approximate global stationary 



5 



2.2 Post linear or second order approximation (exterior solu- 
tion) 

(2) (2) 

By constructing N a ^ and H a from (14) and solving equations (13) for vacuum 
asuming the stationary and axial symmetries, the following expression for the metric 
to second order is obtained (the calculations has been performed using computer 
algebra) 



y [2] 



+ 



+ 



+ 



+ 



i+ 2Ml 



2M, 



(1)2 



2M ( $ 1) 4 1) 



(1) R (1)' 



M^'B. 



2M. 



[2] 



3r 4 
(i) A W 



AMk A 



1 + 



2M, 



[2] 



+ 



AM, 



(1)2 



3r 4 

(1) 4 (1) 



r „ 2 Jf 1 v 
T2 H ^— Zi 



2J. 



[2] 



3r 2 



3^ 



(i) n (i)' 



M^'B, 



Dn 




(17) 



where we used the notation X^ = X^ + X^ in order to take into account the ho- 
mogeneous part of the solution of (13), which is formally identical to (14). Moreover 
the terms involving 1/r 6 or higher, as well as the terms containing the products of 
J moments by another moment, has been omitted as they are not required for the 
final problem. 



3 Interior solution 

3.1 The energy— momentum tensor 

Let us describe now the stellar model. The energy-momentum tensor describing 
the source is a perfect fluid of constant density it , with pressure p depending only 
on r and 9, i.e. 
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T a[} = (i,u a up + p(g a p + u a u p ) , n = const , p(r, 9) (18) 

where u a is the unit 4-velocity of the fluid. Further, we shall assume that no con- 
vection takes place; so, u a is a linear combination of the two Killing vectors we are 
imposing 

u a =^(C+u( a ) , C = S? , C = 5« (19) 

Finally, we assume rigid rotation (Born), that is, u> constant. Consequently, the 
normalization factor tp has the following expression. 

yj = 1 = (20) 

y -(.9oo + 2ujg 03 r sin 9 + w 2 .g 33 r 2 sin 2 9) 

On the other hand, the Euler equations can be written in this case as 

V Q T Q/3 = -^- = d\ogip (21) 

fi + p 

which for constant density implies 

*="(IH <22) 

where ip s is the value of ip on the surface of the fluid, p = 0. The equation defining 
this surface will therefore reads [2]: 

£ : ip(r, 9) = const = ip s (23) 

Now, in order to handle the energy-momentum tensor in the postminkowskian 
approximation context, it is necessary to asign an order of magnitude to the constants 
appearing in this tensor: /i, lu and First of all we will assume obviously that 
the density is a first order quantity (i.e. linear). Furthcmorc, for reasons which will 
become clear in the last section we will assume that the square of angular velocity is 
also of first order. So, if we perform the following substitutions 

ffoo -> -1 + h(xt , 303 -> ^03 > 333 = fd + fa - h 

^ - i + 



we finally get the following tensor T to second order (in the cylindrical basis) 
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7-p] 



2 

+ f 



1 + - 3^ 1} + -yr 1 sin 2 6 



f 00 



fi ur sin 6* t' 



03 



sin 



1.(1) i 2 2-2/1 

/i ' + -u r sin 6* 



t° + f/4V + f 



2 2-2 

cj r sin 



t 2 (25) 



It should be remarked that terms containing has been neglected even though they 
are postlincar in character. This is in agreement with the remark after expression 
(17). 

3.2 Linear approximation (interior solution) 

Imposing regularity at the origin together with the usual symmetries, the general 
solution of equations (12) for the linear part of the energy-momentun tensor is 



where: 



n=0 



71=1 

oc 



n=0 



n=0 



n=0 



(26) 



E„ =E„ + -(2n+l)D„- -(2n+l)H° 

F„ - -(2n + 1)E„ + F„ - ^n(2n + 1)D„ + ^(2n + l)(2n + 3)H° 



and where hp 1 ^ is any particular solution. The constants ran 1 , jn', an' y have 
no specific meaning as in the exterior case. However, since in the exterior we have 
neglected the multipoles M4, J5, At, i?6 and following, the same will be done for 
the coefficients of the analogous spherical harmonics in the interior. So, we have the 
following linear approximate interior solution 



(27) 



,(1) „■(!) „(1) „ Jl) 



[1] 
Bint 



-1 + m^J T + m£> r 2 T 2 + r Zi + r 3 Z 3 
1 + ro^l D + r 2 D 2 + E + F + 
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+4 1) r 2 E 2 +6< 1) r 2 F 2 



^r 2 (T + D ) - ^wr 3 Zx (28) 




3.3 Post linear approximation (interior solution) 

The postlinear approximation can be constructed in a similar way as for the 
exterior case. This would leads to a similar expression to (17), however it is much 
longer and will not be written here. Nevertheless we shall write down below the final 
expression after the corresponding approximation has been made and the matching 
with the exterior completed. 

4 The matching 

4.1 Expansion of the constants and matching surface 

This subsection is the key to understand the analysis carried out in the present 
paper. First we take an adimensional parameter A which essentialy reproduces the 
post-minkowskian approximation carried out above. Then we use another adimen- 
sional parameter SI describing a slow rotation approximation [5] [10] [12] whose 
first term has spherical symmetry ("post-spherical approximation"). Specifically 
the parameters we have used are: 

A) A is the ratio between m, a constant with units of mass, later identified as 
the "newtonian mass" of the source, and a length r measuring the radius of the star 
in absence of rotation (i.e. spherical symmetry). 

B) SI is the ratio between the centrifugal and gravitational forces at the surface 
of the star. It is therefore proportional to the rotational velocity of the source and 
measures its deviation from spherical symmetry. This parameter is used to perform, 
at each order of A, a second approximation (slow rotation). 



Notice that the exponents of A for the Z„ coefficients will be half-integers. Never- 
theless, with respect to the post-minkowskian approximation the order of magnitude 
of the corresponding terms will equal the integer part of those exponents. The final 
approximation will consist on imposing A 5 / 2 — > and fi 4 — > . 

The fundamental hypothesis of this work is demonstrated in the following ex- 
pansions for the "exact" constants M n , J n , A n and B n (the formal meaning of the 



< 



A = — : Post — minkowskian approximation ? 
r o V 

ujr = A 1 / 2 ^ , £1 : Slow rotation approximation 




(29) 
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term "exact" in our case is X n = linv^oo X„ ). In general, these are expansions 
in powers of A, each coefficient of which is, in turn, expanded in powers of 0. This 
hypothesis is based on the newtonian analysis of a selfgravitating rotating object, 
which has, as is well-known, the Maclaurin spheroids as possible solutions [9]. 



Mo 
ro 



= A 



M 2 

,3 



1 + m^q 2 + o(Q 4 ) +a 2 + M^- 2, n 2 + o(n 4 ) +o(a 3 ) 



(2,0) , » /f (2,2) r)2 



M 2 (1 ' 2) + 0(ft 2 ) 



+ x 2 n 2 



M 2 (2 ' 2) + o{n 2 )] + 0(A 3 ) 



Ji 



= A 3 / 2 J?- 1 * + j} 1 ' 3 ^ 2 + 0(ft 4 ) +0(A 5 / 2 ) 



| = A 3 / 2 3 \4 h3) + 0(n 2 )} + 0(A 5 / 2 ) 



(30) 



^ = \n 2 

r 3 

-y = A 

r 



4 1,2) +0(fi 2 ) 



+ A 2 fi 2 



4 2,2) + 0(H 2 )1 + (9(A 3 ) 



B 



(1,0) 



+ a 2 



(2,0) 



Bi 22) n 2 + 0(ft 4 ) +0(A 3 ) 



^ = A0 2 
^o 5 



(1,2) 



o(r! 2 



A 2 2 



(2,2) 



+ 0(Q 2 ) +0(A 3 ) 



(31) 



For M and i?2, the first terms are independent of ft because they already appear 
in spherical symmetry [6] [11]. For J\ (angular momentum) we can assume that the 
lowest term is linear in f2, which corresponds to a sphere put into rotation. However, 
M 2 and J 3 represent true deformations of sphericity due to rotation; therefore, we 
assume that the expansions start as indicated. Finally, for A 2 and B 4 , we assume a 
similar structure to M 2 ; this is based on the fact that these constants are associated 
to the spherical harmonics of the same order as the quadrupole moment. 

Regarding the constants m, j, a and b in the interior solution, we take similar 
expansions as above. For shortness they will not be included here. 

In agreement with the "post-spherical approximation" we make the suplemen- 
tary hypothesis that the equation determining the surface £ of the star has the 
structure: 



r = r Q +r 2 Q 2 P 2 (cos 6) + 0(n 4 ) (32) 

so that, for vanishing rotation, it becomes the equation of a sphere of radius r . In 
the following subsections we analize the matching problem. Let us advance here the 
solution we find for this surface, which reads: 



10 



J. Martin, A. Molina and E. Ruiz 



1 . ±^o- />,;.,, sW: . o,:<) ; : , ,/ ^ (33) 



4.2 Results of the matching (exterior solution) 

Imposing the condition that the metric (interior and exterior) is continuous with 
continuous first derivatives on the surface of the fluid, the following expansions for 
the constants in the exterior solution arc found: 



^ = A + A 2 (3 + ffi 2 ) +0(A 3 ,r! 4 ) 



M2 



1 



-XQ 2 



69 
70 



x 2 n 2 + o(A 3 ,o 4 ) 



Jl 

r 2 

■h 



a : » m ( 1 ■ !<>•' )-o ■ Mn 



i 



A 3 / 2 tt 3 + 0(A 5 / 2 ,r! 5 ) 



(34) 



A 2 
Bo 



(35) 



In the M expansion, the different contributions (newtonian mass, binding energy, 
rotational energy) to the total mass of the body can be observed. 

Plugging this expressions into (17) the following metric tensor of the exterior 
region is obtained 



[2] 
Scxt 





( 2fl 2 \ 




2A 2 ~ 




_l + ^ + 2A 2 








T + ( 


V 


v 3 + - 


; 


T] 2 



A 69A 2 



+ 



2A 5 



+2A3,2s! (f4)? z ' 



2A 3 / 2 » 3 
7n 4 



2A , / 


2^! 2 ^ 




4A 2 ~ 




H h 2A I 3 H 








D + ( 


v V 






" V. 



2r>2 



A 2 ^ 
7r/ 4 

+ A 2 



E 2 + A 2 



1 4 
V + 35 



Q 2 



- <)2 ( 1 



V 14r? 4 63r) 5 



rj 3 35t7 3 n 4 



A 69A 2 32A 2 
~ ~ 35^ ~ 21rf 



Sr T 2 



2r! 2 ^ 



r> 2 d 2 



3 / r? 3 2I77 4 



(36) 
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We should remark that, as mentioned in the abstract, the Kerr metric does not belong 
to this approximate solution. Kerr's metric can be obtained from (17) by making the 
substitutions Mq = m, M 2 — —ma 2 , J\ = ma and J3 = -|mo 3 (here m and a stand 
for the standard parameters of the Kerr metric) . It can be easily checked that these 
values can not be obtained by any choice of the free parameters A, O and tq. One 
might think that this need not to be true for another fluid, for instance a polytropc 
or differential rotation. 

4.3 Results of the matching (interior solution) 

Similarly, the matching conditions produce the following expressions for the 
constants of the interior solution 



mo 



3A + |a 2 (21 + 20 2 ) + 0(A 3 , 4 ) f 3l r = A 3 / 2 (2 + ^O 2 ^) + 0(A 5 / 2 , fi 5 ) 



m 2 r 2 = -AO 2 - t^A 2 2 + <3(A 3 , O 4 ) 



j 3r 3 = _^A 3 / 2 3 + 0(A 5 /2 j f]5) 



(37) 



a = ^A 2 (21 + 20 2 ) + 0(A 3 , O 4 ) { b = 0(A 3 , O 4 ) 
a 2 rl = 



43. 
15 



A 2 2 + 0(A 3 ,O 4 ) 



86 (38) 
fo 2 r 2 = -^A 2 O 2 + O(A 3 ,O 4 ) 



so that the metric in the interior region once matched to the exterior reads: 



[2] . 
Bint " 



-1 + 3A - \r) 2 + \x 2 (| + O 2 ) + A 2 (-1 + O 2 ) n 2 + ^A 2 (\ - \n 2 ) , 4 



+ (-a-^a 2 + ^av)oVt 2 + a 3 / 2 



n 2^ 2 6 2 
2 + -0 2 --, 2 



fi?;Zi - -A 3 / 2 OV Z 3 



+ 



1 + 3A - V + ±A> § + Irf) - + Jrf) + J* I - ) ,< 



D r 



- ^AV + |aV) O 2 D 2 + 1a 2 2 (-4ry 2 + 3r? 4 ) E 2 



+ A 



l(-l + ^), 2 + A ( l-40 2 ), 4 



F 2 + I ^AWF 4 



(39) 
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On the other hand the expression for the pressure is the following: 
A 2 



P : 



8irr 2 



(3-2n 2 )(l-n A )~5ft 2 r 1 2 P 2 (cos9) + 0(A J ,fT) (40) 
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